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Abstract 

We have developed an efficient algorithm for the subtraction of infrared divergences that arise in 
the evaluation of QED corrections to the anomalous magnetic moment of lepton {g—2). By incor- 
porating this new algorithm, we have extended the automated code-generating system developed 
previously to deal with diagrams without internal lepton loops (called q-type), which produced 
convergent integrals when applied to diagrams that have only ultraviolet-divergent subdiagrams 
of vertex type. The new system produces finite integrals for all q-type diagrams, including those 
that contain self-energy subdiagrams and thus exhibit infrared-divergent behavior. We have thus 
far verified the system for the sixth- and eighth-order cases. We are now evaluating 6354 vertex 
diagrams of q-type that contribute to the tenth-order lepton g — 2. 

PACS numbers: 13.40.Em, 14.60.Cd, 14.70.Bh, ll.15.Bt, 12.20.Ds 
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I. INTRODUCTION 



The magnetic moment anomaly of electron, also called the electron g — 2, has been mea- 
sured with steadily increased precision since its discovery in 1947 The latest experiment 
using Penning trap with cylindrical cavity updates the best previous value ^ of the electron 
^-2 to 3 

ae = 1 159 652 180.85 (76) x 10"^^ [0.76ppb]. (1) 

Here, tte = ^{g — 2) and the numeral in the parenthesis represents uncertainty in the last 
two digits of the value. The uncertainty has been reduced by a factor 5.5 over the previous 
best value. 

Theoretically, the electron g—2 is almost entirely accounted for by the QED corrections, 
and thus it has provided the most stringent test of the validity of QED. The dominant 
contribution comes from the interaction between photons and electrons alone (referred to as 
mass-independent term), and it is given as a function of the fine structure constant a. An 
important by-product of the study of the electron g—2 is that a very precise value of a can 
be obtained from the measurement and the theory of Og, assuming the validity of QED. The 



value of a thus determined is {4] 

a-\ae) = 137.035 999 070 (12)(37)(90) [0.71ppb], (2) 

where uncertainties are due to the term, an educated guess of the term j^l, and the 
experiment ([T]). The value ([2]) is by far the most precise determination of a at present^. 

As is seen from Eq. ([2]) the measurement uncertainty in a{ae) is only a factor 2.5 larger 
than that of the theory which is mostly due to the unknown term. When the experimental 
data improves further, the term will become the largest source of unresolved systematic 
errors. This is why an explicit evaluation of term is urgently needed. 

The mass-independent contribution to the term involves 12672 Feynman diagrams, 
which can be divided into 32 gauge-invariant sets. We classify them into 6 super sets (Set 
I- VI) according to their structures (see Figs. I-VI of Ref. |8|). The largest and the most 
difficult is the Set V, which consists of 6354 diagrams that have no lepton loops (we call 
them q-type) and form a single gauge- invariant set. The difficulty of the Set V stems from 



^ Note that the value given in Eq. ([2j includes small changes in the tenth-order estimate [5| and in the 
hadronic light-by- light contribution^, @] , and thus deviates from a^^ given in Ref. 



the fact that many of them, besides being gigantic, have very large number of ultraviolet 
(UV) and infrared (IR) divergences. 

In the previous approach employed in the calculation of sixth-order and eighth-order 
diagrams the subtraction terms of UV divergence were constructed by a procedure called 
il'-operation [ol, which is based on a simple power counting at a UV-divergent singularity. 
In order to deal with UV divergences, let us regularize each relevant photon propagator by 
introducing the Feynman cut-off 



1 11/-^' dL 



(3) 



where A and A are UV-cutoff and IR-cutoff, respectively. Suppose the diagram Q has a 
UV-divergent vertex subdiagram S. Let Mg be the magnetic moment contribution of Q. 
Then the operation of on Mg creates an integral K^Mg which has the same UV cutoff 
A as Mg itself and yet factorizable into the product shown on the right-hand side: 

¥.sMg = (4) 

where L^^ is the leading UV-divergent part (referred to as the most-contracted term) of the 
vertex renormalization constant Lg, and Mg/g is the magnetic moment projection for the 
reduced diagram that is obtained from Q by shrinking iS to a point. It is important to note 
that the factorization in Eq.(j4]) does not work unless both sides are well-defined integrals 
(made finite by the Feynman cut-off or some other regularization). Throughout this paper 
let us assume that all UV-divergent integrals are regularized by the Feynman cutoff. Of 
course the Feynman cutoff is not needed for convergent integrals, and the limit A ^ oo 
must be taken after the renormalization is carried out. 

In Ref. jl^ it was shown that this scheme of UV-subtraction can be incorporated into an 
automating algorithm without change. In the case of Set V diagrams of tenth-order we were 
able to obtain FORTRAN codes for 2232 vertex diagrams containing only vertex subdia- 
grams (represented by 135 Ward-Takahashi-summed diagrams [see Eq. (E])]) out of the total 
of 6354 vertex diagrams (represented by 389 Ward-Takahashi-summed diagrams). A pre- 



shows that UV divergences 



hminary numerical evaluation of these integrals by VEGAS 
are removed completely. 

In case where 5 is a self-energy subdiagram inserted between consecutive lines i and j of 
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Q, the action of on Mg yields a somewhat more comphcated factorization 

KsMg = 6m^/ Mg,si^*) + BJ" Mg/ysj]. (5) 

where 5m^^ and Bg^ are the leading UV divergent parts of the mass renormalization 
constant 5ms and the wave-function renormalization constant Bs-, respectively. The * in 
Mgis{i*) indicates that it has a two-point vertex between the lines i and j. Q /[S,j] denotes 
a diagram obtained by shrinking both S and line j to points. 

This method works as far as subtraction of UV divergence is concerned. However, it 
complicates the handling of IR divergence because of the fact that the unrenormalized Mg 
as well as its self-mass counter term 5msMgis[i*) have a linear or worse IR divergence. In 
the calculation of the eighth-order case we encountered such linear IR divergences in two 



This problem 



diagrams Mig and Mig as well as in some renormalization constants 
has been handled by subtracting the linear IR divergence by an ad hoc subtraction term. 

However, such an ad hoc approach will become very complicated in the tenth-order case. 
Furthermore, it will present a severe obstacle to automation. Thus, we have developed 
an alternative approach to subtraction of the linear (and worse) IR divergence which is 
entirely systematic and fits well in the scheme of automated code generation. It generates 
FORTRAN codes of the renormalized and finite amplitudes for all diagrams of the Set V 
very rapidly, which are ready for numerical integration. 

The aim of this paper is to present the new approach for constructing the IR subtraction 
term in an automated manner. We also report its implementation. 

The organization of the paper is as follows. In Sec. [TTl we briefly summarize our for- 
mulation of the numerical evaluation of lepton g — 2 and the scheme of subtractive UV 
renormalization. In Sec. IIIII we present a new scheme for identifying IR divergences that 
may appear in a complex manner. The construction of subtraction integrals for those diver- 
gences is described in Sec. [IVl In Sec. |V]we show the concrete procedure of code generation 
and its implementation. Sec. |Vl]is devoted to conclusion and discussion. An example of the 
identification of IR subtraction terms is given in Appendix |Al 



^ Mi6 and Afig share the structure which is obtained by inserting two self-energy subdiagrams, one is of 
the second order and the other is of the fourth order, into the second-order magnetic moment term. 
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II. UNRENORMALIZED AMPLITUDE AND UV DIVERGENCES 



In this section we briefly summarize our formulation for evaluating QED contribution 
to the anomalous magnetic moment of leptons by numerical means jo], [l^ • It involves 
the construction of amplitudes for Feynman diagrams and the renormalization of ultraviolet 
divergences. We note that the substructure called forests plays a crucial role in organizing 
the UV renormalization and also accounts for our subtraction scheme of infrared divergences 
discussed in the later sections. 



A. Anomalous Magnetic Moment of Lepton 

The anomalous magnetic moment is given by the static limit of the magnetic form 
factor that is related to the proper vertex part V^. We evaluate the QED contribution to 
in the framework of perturbation theory by the series expansion with respect to a/7r, where 
a is the flne structure constant. 

In our formulation we employ a relation derived from the Ward-Takahashi identity 



A^(p, q) ~ -q 



dAt,{p,q) 



'-^ (6) 



dqy 

between the self-energy part and the sum of vertex parts A'^(p, q) obtained by inserting 
an external vertex in the lepton lines of S in all possible ways. Here, the momentum of the 
incoming lepton p —\q and that of the outgoing lepton p+\q are on the mass shell so that 
p and q satisfy p^ = m? — \q^ and p - g = 0. Eq. ([6]) enables us to turn the evaluation of a set 
of vertex diagrams into the evaluation of a single self-energy-like diagram. This approach 
reduces the number of independent diagrams substantially. 

The magnetic moment anomaly of a 2nth-order diagram Q is given by an integral over 
loop momenta of a product of vertices and propagators of lepton and photon. It is turned 
into a parametric integral over Feynman parameters z-i by the Feynman integration formula. 
Carrying out the momentum integration analytically, we can express the resulting amplitude 
in a concise form: 



Mf) = (--) {n-l)\l{dz) 



, (7) 
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where {dz)g = Yl d^i^i^^'^i ^i)- "^^^ factor (a/vr)" is omitted for simplicity. It is implicitly 
assumed that the Feynman cutoff ([3]) is introduced whenever it is necessary. 

The quantities C^, E^, N^., and are polynomials of symbols called building blocks 



Bij, Ai, and Cij 13j. The symbols Bij and U are homogeneous polynomials of Feynman 
parameters, related to the flow of loop momenta on the diagram. The symbol is called 
scalar current that is associated with the flow of the external momentum p. It is constructed 
from Bij, U, and {zi}. The symbol Cij is also constructed from Bij, U, and {zi}. The symbol 
V in the denominator is a function defined by 

V = J2^^-G, G = J2^^^i^ (8) 

i i 

where the summation is over the lepton lines only and the rest mass of the lepton is set to 
m = 1 for simplicity. 



B. Subtraction of UV Divergences 

The amplitude constructed above is UV-divergent except for the case n = 1. To achieve 
pointwise cancellation of these divergences we adopt here the subtractive renormalization. 

The UV divergence arises when some of the loop momenta go to infinity. In the Feynman 
parametric space it corresponds to a particular regime in which the sum of relevant Feynman 
parameters tends to zero. We prepare the subtraction term in the form of an integral such 
that it cancels the singularity of the original integrand at a singular point of the integration 
domain. These subtraction integrals are constructed from the original integrand by the K- 
operation. As was noted already they factorize exactly into a product (or a sum of products) 
of lower-order magnetic projection and a UV-divergent part of the renormalization constant. 
(See Eq. (H and Eq. p.) 

In the prescription described above, we only subtracted away a part of renormalization 
constants. This is because the integral becomes highly intractable if those renormalization 
constants are treated as a whole. To carry out the standard on-shell renormalization, we 
have to take account of the difference, e.g. Ls — in Eq. (jl]). Similarly for the self- mass 
term and the wave-function renormalization term. We adopt two-step renormalization, and 
adjust the differences afterward by collecting all contributions over the diagrams. We call it 
the residual renormalization step. 
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C. Forests 

A diagram may have complicated divergence structure due to nested singularities. The 
whole structure of UV divergences is managed by Zimmermann's forest formula. A forest 
is a set of divergent subdiagrams and each forest corresponds to a particular emergence of 
divergence. 

First we introduce an inclusion relation of subdiagrams as follows. If two subdiagrams 
share neither vertex nor line, they are called disjoint. If a subdiagram S is included in the 
other subdiagram S' in the sense that all vertices and lines of S are also the elements of 
S', S is included in S'. When S and S' have common vertices and/or lines but one is not 
included in the other, they are overlapping. In this case there are some vertices and lines of 
S that are not the elements of S'. 

A forest / is defined as a set of subdiagrams of which any two elements arc not overlapping. 
It is called normal forest when the diagram Q itself is not an element of the forest. The whole 
set of (normal) forests is found by generating all possible combinations of subdiagrams and 
disregarding any combinations in which a pair of elements are overlapping. 

Noting that the if-operation can be defined for each <S, we can define the UV-finite 
amplitude Mg by the formula 



fed Sef 

where the summation is taken over the set ^ of the normal forests of the diagram Q, and 
the product means the successive application of iC-operations. 

III. IR DIVERGENCES 

The magnetic form factor is free from UV- and IR-divergences once it is fully renormal- 
ized. However, individual diagrams suffer from IR divergences which cancel out only after 
all diagrams are combined. 

The root cause of IR divergence is the vanishing of the denominator of the photon prop- 
agator l/k"^ in the limit k ^ 0. This is, however, not the sufficient condition since it gives 
a finite result on integration over the 4-dimensional momentum k. In order that it becomes 
divergent, it must be enhanced by vanishing of the denominators of at least two lepton 




(9) 
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propagators due to some kinematical constraints. Typically, this happens when the momen- 
tum of each of these lepton propagators is constrained by sharing a three-point vertex with 
the soft photon and an external on-shell lepton line. When the external momentum p is 
constrained by the on-shell condition p'^ = m?^ the lepton propagator in question behaves as 

\ = (10) 

{p + ky -m"^ 2p - k + k"^ 2p ■ k 

for k —>■ 0. These lepton propagators will be called "enhancers". The logarithmic IR 
divergence takes place when the /c-integration is carried out and the soft photon singularity 
is assisted by two enhancers. When the vertex Feynman diagram G{k) in question has a 
self-energy subdiagram, we find three enhancers due to the kinematical constraint of the two- 
point vertex so that we find the IR divergence to be linear. The IR divergence becomes even 
severer when the diagram G{k) has more than one self-energy subdiagram which effectively 
bring in a number of two-point vertices. 

To handle the IR divergences, we adopt again a subtractive approach in which an integral 
of IR subtraction terms is constructed in such a way that it cancels out the IR-divergence 
of the integral Mg'^^ of Eq. (j7]) point-by-point in the Feynman parameter space. 

First we briefly summarize a scheme called /-operation which has been developed and 
employed for the calculation of sixth- and eighth-order calculations. 

In the Feynman-parametric representation, the IR divergence is caused by vanishing of 
the denominator function V of Eq.® in the corner of the integration domain characterized 
by 

Zi = 0{6) if i is an lepton line in TZ, 

Zi = 0{1) if i is a photon line in TZ, 

Zi = 0{e), er^6^ iiieS, (11) 

where TZ = Q/S. In this limit the denominator V vanishes as 0{5'^). (The last condition is 
actually an artifact of the condition Zi = 1, which can be readily lifted.) 

If two lepton propagators participate in the enhancement, we obtain a logarithmic IR 
divergence. In this case we can construct an IR subtraction term by a simple power counting 
rule and an /-operation similar to the /^-operation of the UV divergent case. For the 
subdiagram TZ = Q/S the /-operation I-ji is defined as follows: 

(a) In the limit (fTTj) keep only terms with lowest power of e and 6 in U, Bij, Ai. 



(b) Make the following replacements: 



U^UsUn, V^Vs + Vn, F ^ Fo[Ln]Fs, (12) 

where Fo[L7^] is the non-contracting term of the vertex renormalization constant defined 
on TZ, and Fg is the product of 7 matrices and Df operators for the diagram S. (See 
Ref. jsl for definitions.) 

This procedure creates an integral defined on the parametric space of Mg. By construction 
it factorizes as 

Lnik) [Fo] Ms + Mn* [I] ASrhs . (13) 

For precise definitions of quantities quoted above see Ref. 9|]. 

In the following, we see that for the lepton g—2 of the q-type diagrams, there are two kinds 
of sources of enhancement. These two types of divergences may occur simultaneously in a 
complicated manner. We will introduce a new scheme modeled on the previous approach 
for identifying the emergence of the divergences. In this scheme the identification relies 
on examining the shape of the diagram that allows diagrammatic treatment systematically 
without the need for close examination of the integrand, and it is also suited for automated 
treatment. 



A. IR Divergence Caused by Residual Self-mass 

One type of IR divergence appears as a consequence of our particular treatment of the UV 
divergences by means of K-operation. Suppose a diagram Q has a self-energy subdiagram S. 
As is readily seen from the analysis of Feynman diagrams, this divergence is not the source 
of real problem since it must be canceled exactly by the mass-renormalization counterterm 
6ms Mg/s (i*), where Srris is the (UV-divergent) self-mass associated with the subdiagram 
S defined on the mass shell. The reduced magnetic moment amplitude Mg/5(j*) is the one 
that has a linear IR divergence. As a consequence 

Mg -6msMg/s{i^) (14) 

is free from linear IR divergence. Although this cancellation is analytically valid, however, 
it is not a pointwise cancellation in the domain of Mg. Our problem is thus to translate the 
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second term into a form which is defined in the same domain as that of Mg and cancels the 
IR divergence of Mg point-by-point. 

Now, as was noted in Eq. the i^-operation for the subdiagram S acting on Mg creates 

KsMg = 5mUV Mg/s + 5^ Mg/[s,j] . (15) 

If we find an integral what causes pointwise cancellation of the linear IR divergence in the 
domain of Mg and also produces the factorization as 

SmsMg/s(i*) (16) 

where 

Sms = Smg — Srn^^ , (17) 
then from Eqs. f|T5l) and fflGl) we would have 

KsMg + 5^sMg,s ii*) = 5ms Mg/s (i*) + B^^ Mg/[sj] ■ (18) 

If we schematically introduce an operator that produces the integral of Eq. f|T6l) as 

RsMg = 5^sMg/sii*), (19) 

Eq. (|T8l) would then be written as 

{Ks + Rs)Mg = 5ms Mg/s (i*) + Bf^ Mg/ys,j] • (20) 

It turns out that it is not difficult to construct such an integral. Furthermore, it can be 
readily incorporated in our automation algorithm. We call this subtraction scheme as the 
residual self-mass subtraction, or "i?-subtraction" operation. 

B. Modified /-subtraction operation 

After the linear IR divergences are disposed by the iT-operation and /^-subtraction op- 
eration we are still left with logarithmic IR divergences. To treat these divergences let us 
consider a vertex diagram Q{k) which has a subdiagram S{k). Here k refers to an external 
photon vertex attached to a lepton line Ck of S. The reduced diagram IZ = Q{k)/S{k) is 
connected to S{k;) by lepton lines and £j. (See Fig. ^l(a)[ ) 
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This diagram exhibits an IR-divergent behavior when the momenta of (all or some) 
photons in IZ go to zero, accompanied by the enhancers Ci and Ij. The substructure S{k) to 
which these enhancers are attached can be considered as a magnetic moment of the lower 
order. Thus the amplitude in this limit becomes that of the diagram TZ{k) obtained by 
replacing S{k) by a vertex that is weighted by Ms(k), as shown in Fig. [ ^b) 7^(/c) develops 
a logarithmic IR divergence as is easily verified by power counting. 

Since we are dealing with the Ward-Takahashi-summed diagram defined by Eq. ([6]), we 
have to consider the sum of contributions of vertex diagrams S{k) which are obtained by 
inserting an external vertex k to the self-energy-like diagram S in every possible way. In 
the calculation of the diagram the IR singularity associated with the substructure S 
has a form that consists of contributions from the vertex diagram TZ{k) and the magnetic 
projection of the diagram S. The IR singularity is contained in the vertex renormalization 
constant L-jK^^y 

For the explicit expression of the IR subtraction term that cancels the singularity in the 
above, we consider the term 

Lg/sik)Ms, (21) 

where L is the residual part of the vertex renormalization constant 

Lg/sik) = Lgis{k) - Lgjs(k) ■ (22) 

We construct an integral that corresponds to Eq. ( l2Ti) in the domain of Mg. We call this 
subtraction scheme as the "/-subtraction" operation. 
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Note that the subtraction scheme for IR divergence adopted here is different from that 
of Ref. [QJ. We now choose as the IR subtraction term the part of Lqisq^^ which excludes 
the overall UV divergent part but includes a finite part /S.Lgis{k) left out in the previous 
method. This simplifies the handling of the IR problem considerably. 

The magnetic moment part Mg may have UV and IR divergences, which also have to be 
subtracted. The UV sub divergences of Ms will be treated in a similar manner to the UV 
sub divergences of L. The IR divergence that is related to the nested singularity is discussed 
in the next subsection. 



C. Nested Singularity 

We have introduced the operations for identifying the subtraction terms called /- and 
i?-subtraction operations corresponding to two types of divergences. There may appear 
more than one source of IR divergences that lead to complicated divergence structure. They 
are treated by combinations of /-/i?-subtraction operations conducted by annotated forests 
described later. 

By generalizing Eq. (|T9|) the i?-subtraction operation for a self-energy subdiagram S 
is defined by 

^sMg = 5m'lMg,sii^): (23) 
where the operator R5 acts on the UV-finite amplitude M_g of the self-energy-like diagram 
Q and generates a product of a residual part of mass renormalization constant 5m^ defined 
by Eq. (1211) and a UV-finite lower-order amplitude Mg/5(j*)- The reduced diagram Q /S {%*) 
is obtained from Q by shrinking 5 to a point. The term is defined as a residual part of 
the mass renormalization constant 5ms defined on the mass shell by subtracting the leading 
UV divergence (referred to as the most-contracted term) 5m}g^ and the sub divergences as 

5m^ = 5ms - 5mJ^ + 11 ^-^s') ■ (24) 

/ s'ef 

5ms = 5ms — 5mb^ may also have extra UV divergences due to its substructures that are 
subtracted by the iC-operations as above. Here, the summation is taken over the normal 
forests of the subdiagram S. 

The /-subtraction operation I5 for a self-energy subdiagram S is defined by 

lsM^^Ll,sik)M^. (25) 
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where the operator I5 acts on the UV-finite amphtude M g of the self-energy-hke diagram Q 
and generates a product of a residual part of vertex renormahzation constant Lgy^j-^-j defined 
by Eq. ( l26l) and a UV-finite lower-order amplitude M^. Here, the diagram lZ{k) = Q /S{k) 
is given from Q by replacing S by an external photon vertex labelled by k. The term L^^^-^ is 
defined as a residual part of the vertex renormahzation constant L-jk^^-^ defined on the mass 
shell by subtracting the leading UV divergence (referred to as the most-contracted term) 
L^y^-j and the sub divergences as 

Ln(k) = L-R-ik) - L^k) + X] n (^^5') L-TK^k) ■ (26) 

/ 5'G/ 

L-jK^k) = -^7^(A;) ^ Ln(k) have extra UV divergences due to its substructures that 

have been subtracted by the iT-operations. The summation in Eq. fl26|) is taken over the 
normal forests of the diagram TZ{k). 

For a diagram Q containing a single self-energy subdiagram S, the associated IR diver- 
gences are treated by those two types of operations, and the IR-finite amplitude is thus 
given by 

AMg = Kg- IsMg - ^sKg . (27) 

When the diagram has more than one such self-energy subdiagram, the IR divergences due 
to all those subdiagrams have to be subtracted away. The finite amplitude free from both 
IR and UV divergences are provided schematically by 

AMg = J](l-l5-K5)Mg, (28) 

s 

where the product is taken over all self-energy subdiagrams of the diagram Q. 

By expanding the product in Eq. ( l28l) . we are lead to a forest-like structure which is 
analogous to the renormahzation of UV divergences. In this case, a forest consists of only 
self-energy subdiagrams, and each subdiagram is assigned a distinction of i?-subtraction or 
/-subtraction. We call such a forest "annotated forest" . Eq. f l28l) is thus turned into a sum 
over all annotated forests / as 

^(-l5j---(-M5,)---Mg, (29) 
/ 

where iSj, . . . and Sj, . . . are elements of the annotated forest / that are assigned to I- 
subtraction operations and /^-subtraction operations, respectively. Those operators act on 
the amplitude M g successively in an order given below. 
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FIG. 2: Successive operations of two /-subtractions for subdiagrams Si and 1S2 that satisfy Si D 
They yield a product of component terms shown on the right-hand side. 




By construction of /- and i?-operators, there are some annotations that are not accepted. 
To see which annotations are allowed, it is sufficient to examine the cases in which two self- 
energy subdiagrams are concerned. The inclusion relation of two self-energy subdiagrams is 
disjoint or inclusive. Therefore, the possible patterns are exhausted by the cases in which 
both subdiagrams are assigned to /-operators, both subdiagrams are assigned to /^-operators, 
or one subdiagram is assigned to /-operator and the other to /^-operator, for two self-energy 
subdiagrams which are disjoint or inclusive. We elucidate on those cases separately as 
follows. 

First we consider the case which involves two /-subtractions. Recall that the action of 
operator I^^ is given in Eq. (!25|) . It is found that an extra operation of operator I^^ for the 
subdiagram ^2 must have its domain in M^^. Therefore the subdiagrams have to respect 
the relation Si D 52. Otherwise, the operation of I^^ turns to zero. It leads to a rule that 
subdiagrams for /-operators must satisfy that one subdiagram is included in the other, and a 
pair of /-operators for disjoint subdiagrams 5i niS2 = is not allowed in the annotation. The 
ordering of /-subtraction operators is given by Is2^Si for Si D ^2. (The outer subdiagrams 
are operated first.) 

h,hMg = L^/Siik)Ll/s,(k')Ms, for Si D ^2. (See Fig. El) (30) 

The /^-operator for a subdiagram Si acts as Eq. (l23l) . If an extra operator of is 
applied, ^2 and the reduced diagram Q/Si (i*) must share some parts. Thus Si is included 
in 1S2, Si C ^2, or Si and ^2 are disjoint. The ordering of two /^-operators is given by'^ 

^ The operation of this case is actually recognized as K^^/g^. 
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FIG. 3: Successive operations of two iZ-subtractions for subdiagrams <Si and ^2 that satisfy Si C S2- 

g ^ --^QIS'i So/Si Si 



FIG. 4: Successive operations of two iZ-subtractions for subdiagrams Si and 52 when they are 
disjoint, 5i n ^2 = 0. 




for Si C ^2. (The inner subdiagrams are operated first.) 

M&M^iMg = Mg/s, i^*J^l/s, i^l)Srnl for Si C ^2. (See Fig. O) (31) 

For disjoint subdiagrams, the ordering is indifferent, 

^s^^sMg = ^Si^sMg = Mg/(5iU5.) i^lq)S^sMl for 5i n ^2 = 0. (See Fig. H) 

(32) 



Next we consider the case in which there are two operators I^j and M^j- Since I^^ 
yields M^^ and yields Mg,^„ ^-^^ Si and G/S2 must share some parts, i.e. Si D ^2 or 
SiH S2 = 0. For the case Si D ^2, the successive operations yield 

^sMKg = Rs, {L^/s, M5J = L^/s, M^xis, iq) ^< , (See Fig. El) (33) 
For the disjoint case iSi fl 1S2 = 0, it is given as 

hi^sMg = iq) 5<) = L^/s,/sAk) ^s, ^< , (See Fig. 0) (34) 

The case with Si C ^2 is prohibited. 

To summarize, an IR-finite amplitude for a self-energy-like diagram Q is obtained by 
subtracting away all IR divergences associated with its self-energy subdiagrams. They can 
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FIG. 5: Successive operations of /-subtraction for subdiagram Si and i?-subtraction for subdiagram 
«S2 when they satisfy Si D 62- 




FIG. 6: Successive operations of /-subtraction for subdiagram Si and /^-subtraction for subdiagram 
^2 when they are disjoint, <Si n (S2 = 0. 




be identified by applying two types of subtraction operations to the amplitude, called I- 
subtraction operation and /^-subtraction operation (M5), for a set of annotated forests 
similar to Zimmermann's forests for UV divergences. An annotated forest consists of self- 
energy subdiagrams, to each of which a distinction of /-subtraction or /^-subtraction opera- 
tion is assigned. By construction, the annotation has to respect the assignment rules: 

• subdiagrams assigned to /-subtraction operations are included one another, 

• a subdiagram S' assigned to /Z-subtraction operation does not include a subdiagram 
S assigned to /-subtraction operation. 

The IR-finite amplitude is thus given by 

AMg -M.g + J2 i-W ■ ■ ■ i-^s,) ■■■Kg, (35) 

in which the summation is taken over all annotated forest / = {Si, . . . ,Sj, . . .}. (Recall that 
M r, is UV-finite.) The allowed combinations of the operators, including their ordering, are 
listed as follows: 
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11 



111 



IV 



Is, ■ for Si C Sj, 

Rs, ■ Rs, for Si D S,, 

Rs., ■ Rs^ = M5, ■ for S, n Sj 

Rs, ■ Is, for Si C S„ 



v) for^^n^, =0. 

It is noted that the operators Is and M5 are regarded to act on the UV-finite amphtude 
M g in which UV divergences are subtracted away by if-operations. The /-/i?-subtractions 
are introduced in such a way that they produce IR subtraction terms as products of UV- 
finite quantities. To be consistent with the formal prescription for successive operations, we 
have adopted the definitions of Eqs. (!23|) and (!25|) . However, the identification of the IR 
subtraction terms can be carried out diagrammatically by the form of the diagram alone, 
and we do not have to examine the internal UV structure of the diagram nor the explicit 
expressions of integrands. Thus we introduce the regulations as a formal treatment that the 
/-/i?-operators satisfy I^K^/ = and = for /-(i?-)subtraction operator I5 {Rs) for 

a subdiagram S and a if-operation K5/ for a subdiagram S'. 

IV. CONSTRUCTION OF IR SUBTRACTION INTEGRALS 

In our numerical approach, the subtraction terms are given in the form of the parametric 
integrals that are defined over the same Feynman parameter space as that of the original 
diagram. These terms are prepared so that the IR singularities of the integrand are canceled 
at the same point. To attain this point-by-point subtraction of the singularity associated 
with the IR limit, we construct the IR subtraction term L^^^f^^M^ in such a way that the 
component terms L^^j^^ and expressed by separateparametric integrals are merged into 
a single integral by the Feynman integration formula 9l. Il0|. 

In the previous section we have seen that the IR subtraction terms can be identified 
by annotated forests of the diagram Q. By the successive applications of /-/i?-subtraction 
operations along the forest, the form of each term is obtained as a product of component 
terms, which are the residual part of the vertex renormalization constants and/or that of 
the mass renormalization constants, and the lower-order magnetic moment part M. 
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In the actual construction, we follow the two-step procedure. First, we construct an 
integral that corresponds to a product of M, L = L — L^^ and/or 5m = 6m — 6m^^ . 
The location of component terms in the integration space, i.e. how the reduced diagrams 
of the components are embedded in the original diagram, is crucial to realize the point- 
wise subtraction. Here we introduce a tree representation of the annotated forest, called 
as an annotated tree, which provides a useful tool for identification of Feynman parameter 
assignment of the component integral. 

In the next step we subtract away the UV sub divergences of those terms by applying 
if-operations that are restricted onto the respective reduced diagrams. For this purpose, 
we introduce an embedding tree of the annotated tree, which corresponds to a particular 
combination of UV sub divergences of component terms. 



A. Subtractive Integral in Feynman Parameter Space 

The IR subtraction term is given as a product of component terms. Each component 
term is expressed as a parametric integral over the reduced diagram. Consider a case in 
which the IR subtraction term consists of two component terms of reduced diagrams S and 
71 = Q/S. Let xi,X2, ...yXng and yi,y2, ■■■illnu be the Feynman parameters of the diagram 
S and 7^, and let xs = Tli^i and yu = J^Viy respectively. Then the IR subtraction term 
takes the form 

W.^x|(dy),^, (36) 

where we denote the numerators of the part of integrands with the particular powers a and 
P of V^-function as g[S] for the component of the subdiagram S, and g[7l] for the component 
of the reduced diagram 71, respectively. 
By inserting the identity 

in Eq. (136|) and making use of the Feynman integral formula 

v(k) rd] „k-ifi-i 
J^A^^r(,u)l.ism^-s-,)^-^^^, (38) 

the product of integrals in Eq. ( !36l) is turned into a single integral over Q of the form 

g[S]g[7Z] 

{Vs + VnY^P ' ^^^^ 
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FIG. 7: An example of nested subdiagrams of a forest and the tree representation. 




where 

Zi = XsXi for i = 1, ... ,ns, 

Zi = ynVi ioi i = ns + 1, . . . , ng. (40) 

By this assignment the singularities in the Feynman parameter space of both the original 
integrand and that of the subtraction term cancel each other. This mechanism of cancellation 
is crucial for our subtraction procedure. 



B. Tree Representation of Annotated Forest 

A forest is represented as a tree form that expresses the inclusion relation of the subdi- 
agrams in the forest. (See Fig. [3) We assign each subdiagram to a node. If a subdiagram 
is included in another subdiagram S, it is expressed as a child node of the node assigned to 
S. We consider the diagram Q itself as the root node of the tree. For later convenience, we 
denote the subdiagrams that contain another subdiagram S as (direct) ancestors of S, and 
the subdiagrams that are included in S as descendants of S. 

For the annotated forest, a distinct operation of /-subtraction or /^-subtraction is assigned 
to the node. We call such a tree as annotated tree hereafter. Each node is then translated 
into the component term of the reduced diagram. 

For a simplest example in which Q has a single subdiagram S, the tree form of the forest 
is shown as follows: 

g s 
© — ► o 
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When we consider the I5 operation, the result is given as a product of a vertex renormahza- 

tion constant -Z^^^^) for the part of the diagram TZ{k) = Q /S{k) and a magnetic moment part 

M_g for the subdiagram S. Those parts of subdiagrams are related to the nodes in the above 

tree labelled by Q and iS, respectively. Thus we represent the assignment of components of 

the subtraction term graphically as below. 

Q S 
© O 

When the i?-subtraction is considered, the assignment is represented in a similar manner 
as follows. 

g s 
@ — > o 

Here, the subdiagram S referred by the right node is related to the residual self-mass term 
Srn^, and the reduced diagram Q/S (i*) that corresponds to the left node is assigned to the 
magnetic moment part M_g/s{i*)- 

This representation can be extended to more general cases. For an annotated forest that 
corresponds to a nested singularity, the successive operations of /-/i?-subtractions in an 
appropriate order as described in Sec. IIII CI is interpreted to extend the tree by following 
the respective process shown above. Then we obtain a tree representation of the annotated 
forest in which the nodes are assigned distinct types of M-, L^-, or 5m^-nodes. 

The nodes of the tree are related to the reduced diagrams, and so the annotated tree has 
a direct interpretation to the IR subtraction term in the form of a product of the component 
terms for their respective reduced diagrams. Since the tree expresses the nested structure 
of subdiagrams in the forest, the relevant set of Feynman parameters for the component 
term, i.e. how the reduced subdiagram for the component term is embedded in the original 
diagram, can be easily read off from the tree. This feature is crucial in constructing the IR 
subtraction integral so that the point-wise subtraction of IR singularities is achieved. 

Thus far, the annotated tree provides a graphical representation of the annotated forest, 
and it has a one-to-one correspondence. It is readily translated into a symbolic form of the 
associated IR subtraction term, which is also significant for the residual renormalization step. 
The set of IR subtraction terms can be obtained by finding the set of annotated trees that 
have proper assignment of types to the nodes of tree consisted of self-energy subdiagrams. 
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The rules for the annotated forests given in Sec. IIII CI are reflected to those for the 
assignment of types of components in the annotated trees. They are summarized as follows: 

1. There is one and only one node to which the magnetic moment part M is assigned. 

2. The nodes that are assigned to the /-subtraction, L^, are restricted to the ancestor 
nodes of the magnetic moment part. 

3. The nodes that are assigned to the residual self- mass subtraction, 5m^, do not appear 
as the ancestor nodes of the magnetic moment part. 

It turns out that to fulfill these rules the assignment is uniquely determined once a node 
is chosen for the magnetic moment part. We first pick up a node that is assigned to the 
magnetic part M, and then the nodes that lie as ancestors of the M-node are associated 
with the /-subtractions. The remaining nodes are assigned to the /^-subtractions. 

C. Subtraction of UV Subdivergences 

We now proceed to the treatment of UV subdivergences in the components of IR subtrac- 
tion terms. For example, the residual part of vertex renormalization constant is given in 
the form 

^nik) = Ln(k) - L^k) + ^ Yl ("^^O ^^k) ■ (41) 

/e5 S'ef 

We have introduced a term L-jK^k) = /^7e(fe) — -^t^II) construction of the subtraction term 

which is given from L-jK^k) by dropping the most-contracted term L^^^^ . The divergences that 
originate from the substructure have to be subtracted away that are dealt with /C-operations 
for the forests / of TZ{k). In this section we show that these subdivergences can be handled 
in a systematic way and present it as a general scheme. 

The IR subtraction term constructed in the previous section is originally given as a 
product of separate integrals defined on the respective component diagrams. As an example, 
the term L^^^-^M_g is made of the vertex renormalization constant defined for the vertex 
diagram TZ{k) and the magnetic moment part defined for the diagram S. Those components 
have their own divergent subdiagrams. 
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FIG. 8: A UV divergent subdiagram S in the annotated forest {Q,S} 




It is shown that these subdiagrams of separate components can be mapped to subdiagrams 
of the original diagram Q. A subdiagram 5 of a component diagram is itself a subdiagram 
of ^, or 5 is obtained as a reduced diagram by shrinking other component diagrams. 

Thus, it can be recognized graphically in a way that a forest of a component diagram is 
related to a forest of the original diagram. The tree form of the forest of the component 
diagram can be embedded as a subtree of the original annotated forest that corresponds 
to the IR subtraction term by mapping the elements of the subtree as that of the original 
forest. We call such a tree that corresponds to UV subdivergence of certain components the 
embedding tree. For example, consider a case in which the component diagram TZ = Q/S 
has a forest consisting of a subdiagram S as shown in Fig. [HI The tree form of this forest 

n = g/s s 
© — — o 

is embedded in the annotated forest 

g s 
© — > o 

which reduces to the forest of the whole diagram g shown as follows. 

g s 

" o 

It is noted that the component diagrams are disjoint by construction and therefore the 
subdiagrams of different component diagrams are not overlapping. 

The UV sub divergences contained in the product of the components are identified by 
forests of the original diagram g. For the annotated forest / that corresponds to the IR 
subtraction term, consider a forest /' that includes /, in the sense that all subdiagrams of 
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/ are also the elements of /'. Then the subdiagrams in /' — /, i.e. the subdiagrams of 
/' that are not in /, are related to the sub divergences of the component diagram to which 
they belong respectively. The subtraction term of the UV subdivergence is constructed by 
applying iC-operations for these subdiagrams restricted onto the corresponding component 
diagram. For example, acts on the component Q/S in Fig. ([8]). 

Therefore, the UV sub divergences for an IR subtraction term that originate from the 
distinct component terms are synthetically identified by the forests of the original diagram 
Q. However, not all the subdiagrams of Q cause divergences of the IR subtraction term 
related to the annotated forest /. Consider a case in which a vertex subdiagram S' includes 
a self-energy subdiagram S, i.e. S' D S. For the /-subtraction associated with S, the 
subdiagram S is replaced by an external vertex. Thus S' acquires the extra vertex in the 
reduced diagram Q/S{k) and becomes irrelevant for the UV divergence. We have to exclude 
such cases from the identification of sub divergences. 

Another rule must be imposed in relation to the residual self-mass subtraction. First 
consider a nested UV divergences in which self-energy subdiagram S is included in a vertex 
or another self-energy subdiagram S. There would be a UV subtraction term denoted as 
5m^^L^^^.^^ for the vertex subdiagram, or 6rn^^ B^J^ ^.^.^ for the self-energy subdiagram in 
the prescription of i^-operation. However, by construction of i^'-operation it is found that 
those terms that involve the most-contracted part of the vertex renormalization constants 
or wave-function renormalization constants with mass insertions do not actually appear. 
Accordingly, we also have to omit the residual self-mass term 6m^L^y^ or dm^B^Y^ 
respectively. They would be realized as UV sub divergences of the component S/S. 



V. PROCEDURE AND IMPLEMENTATION 



In this section we describe the flow of process to generate the numerical integration code 
for IR subtraction terms of a q-type diagram from its representation. 

We first identify the IR divergent parts based on the forest structure of the diagram and 
provide the set of subtraction terms as annotated forests. Next we identify the UV subdi- 
vergences of these IR subtraction terms. The subtraction terms for these sub divergences are 
constructed by applying i^'-operations to the IR subtraction integrals. 

The subtraction terms of IR divergences and their UV sub divergences thus constructed 



23 



are turned into FORTRAN codes. They are merged to the numerical calculation codes for 
UV finite amplitude generated by the scheme developed in Ref. [lO]. The entire fiow of 
process is depicted in Fig. O The numerical integration codes are readily integrated by 
VEGAS to produce a finite amplitude for the diagram. Some of the steps are common to 
the procedure for UV subtraction described in Ref. jlo| . 



A. Finding Forest Structures 

The types of subdiagrams that are relevant to the IR subtraction term and their UV 
sub divergences are self-energy subdiagrams and vertex subdiagrams. For a q-type diagram a 
subdiagram of these types is represented by a segment of the lepton path. All subdiagrams 
of a diagram are thus found by examining all segments that correspond to the connected 
one-particle irreducible subdiagrams of vertex type or self-energy type. 

The inclusion relation of the subdiagrams are thus mapped to that of the segments. Once 
it is found, the whole set of (normal) forests is found by generating all possible combina- 
tions of subdiagrams and disregarding any combinations in which a pair of elements are 
overlapping. 

For each forest we examine the nest structure of subdiagrams and describe it in a tree 
form whose root node is associated with the whole diagram Q. 



B. Finding IR Subtraction Term 

The IR divergences of a q-type diagram emerge in relation to subdiagrams of the self- 
energy type. To identify the divergence structure, we first pick out the forests that consist 
of only self-energy subdiagrams. Suppose they are expressed in a tree form. Next, we assign 
to each node of the tree the type of component of subtraction term among the choices: the 
residual part of vertex renormalization constant {L^) that is related to /-subtraction, the 
residual self-mass part {6m^) that is related to i?-subtraction, and the magnetic moment 
part (M). The assignment that fulfills the rules in Sec. IIIII can be found by following steps. 

1. Choose a node for the magnetic moment part (denoted by M-node). 

2. The ancestor nodes of M-node are assigned to the L^-part (denoted by L^-node). 



24 



3. Other nodes are assigned to the residual self-mass parts (denoted by 5m -node). 

We call the forest with proper assignment of types to the nodes of its tree representation as 
the annotated forest. Each annotated forest corresponds to a IR subtraction term. 

The expression of the IR subtraction term associated with the annotated forest is then 
given in a form of a product of component terms that are related to the nodes of the tree. 
The component term of a node S is found by following the rules below. 

• An L^-node corresponds to the residual part of the vertex renormalization constant 
-^^(fc) diagram TZ{k). It is obtained from S by replacing the inner subdiagram 
that is assigned to M-node or L^-node by a vertex k, and by shrinking the subdiagrams 
of child nodes to points. 

• A (5m^-node corresponds to the residual part of the mass renormalization constant 
Sm^ for the diagram TZ. It is obtained from S by shrinking the subdiagrams of child 
nodes to points. 

• An M-node corresponds to the magnetic moment part for the diagram TZ obtained 
from S by shrinking the subdiagrams of child nodes to points. 

As an example, consider an annotated tree shown as 




It corresponds to the IR subtraction term 

Lg/Sr/S2ik,){h*) S^Sr ^52/53 fe*) ^^ly (42) 

Since the mass renormalization constant of second order 6m2 satisfies (5m^^ = 6m2, the 
residual self-mass subtraction for the second-order self-energy subdiagram is not necessary. 
We omit the annotated forests that involve such cases. 

An example of a complete set of IR subtraction terms for an eighth-order diagram iden- 
tified in the above procedure is shown in Appendix \M 
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C. Constructing IR Subtraction Integral 



The integrand of the IR subtraction term associated with an annotated forest / is con- 
structed by the following steps. First we identify the position of each component term in the 
original diagram Q, and find the respective sets of Feynman parameters for those compo- 
nents. This step is crucial for our numerical method based on the point-by-point subtraction. 
Then we follow these steps: 

1. Construct the integrand as a product of those of component terms. 

2. Find the building blocks, B^, U, Cij and Ai. 

3. Construct the ^-function in a form that guarantees the analytical factorization of the 
subtraction integral into components. 



1. Integrand 

The numerator of the integrand of the subtraction term is given as a product of those of 
component terms. We first consider the full part of the component integrands, the vertex 
renormalization constant L, the mass renormalization constant 6m, and the magnetic part 
M, before dropping other than the residual part. 

For example, consider an annotated forest consisted of an /-subtraction for the subdia- 
gram S. We construct the subtraction term of the form of product 

Ln(k) Ms. (43) 
The numerator of the integrand (except for the numerical constants) is thus given by 

Tr7^(fc) |p.7"(^ -m)---f,- Tr^jpT^l^ - m) ■ ■ ■ • (^^ - m) ■ ■ ■ 7,, | , (44) 

where the labels TZ{k) and S of the traces indicate that the traces are taken separately 
for the components in TZ{k) and S. Py and P are the projection operators for the vertex 
renormalization constant and the magnetic moment, respectively, and the term refers to 
the external vertex k. 

The expression (14^ is analogous to the integrand of the unrenormalized amplitude of the 
original diagram Q except for these points: 
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• There are multiple traces, each of which corresponds to the respective component 
term. 

• Each part has appropriate projection operator inserted. 

• For the L^-parts associated with the /-subtractions, extra vertices are inserted. 

On the other hand, the pattern of contractions of the internal vertices by the photon prop- 
agators is the same as that of the original diagram Q given by the pairings of the diagram 
representation. Thus we construct the integrand of the subtraction term based on that of 
the unrenormalized amplitude by incorporating modifications listed above. 

We perform the trace operations, and turn the integrand into a Feynman parametric 
form by carrying out the momentum integration analytically using the integration table. 
The result is expressed as a polynomial of the symbols called the building blocks, B^j, U, 
Cij, Ai and V, each of which is given as a function of Feynman parameters. 

We have to remove the leading overall UV divergent parts of the component terms L and 
dm to leave their residual parts. It is achieved by dropping the most-contracted term of the 
component integral which can be identified by simple power counting. 

2. Building Blocks 

The integrand of the IR subtraction term is constructed as a product of component terms, 
each of which is defined on the distinct set of Feynman parameters and given in the form 
of rational function of building blocks that are also defined on the respective component 
diagram. However, we need not prepare these building blocks separately if we recall that 
the building blocks B^j and U factorize exactly in the UV limit. 

In the UV limit associated with a subdiagram S, the building blocks Bij and U factorize 
in the form 




(45) 



for i,j e S, 



(46) 



for i,j e g/S, 



(47) 



g Q/s 
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where the labels Q, S, and Q/S denote that the quantities are defined on the respective 
diagrams. 

Therefore, we are able to use the building blocks Bij and U obtained in the UV limit 
associated with the same forest /. We have to eliminate unwanted contributions by explicitly 
putting the elements of Bij to zero for i G iS and j ^ S to ensure that the integrand originates 
from the separate integrals. 

3. V -function 

In order that the exact factorization of the integrand holds, the l^-function in the de- 
nominator of the integrand is constructed in the form 



where the sum is taken over the component diagrams. 
D. Subtraction Terms for UV Subdivergences 

The integral forms of the subtraction terms for UV subdivergences that originate from the 
components of the IR subtraction term are constructed by the following steps. We assume 
that the IR subtraction term is associated with the annotated forest /. 

1. Identification of Subdivergences 

The set of UV subdivergences to the annotated forest / is given by finding forests /' 
that include /. The emergence of UV subdivergences are associated with the subdiagrams 
in (/' — /), i.e. the subdiagrams of /' that are not the elements of /. By the tree forms 
of the forests /' and / we identify which component of IR subtraction term each of these 
subdiagrams belongs to. We check if these subdiagrams in the component diagrams actually 
cause the divergences, and disregard the cases described in Sec. IIV CI 




(48) 



s 
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2. Integrand 



The integrand of the subtraction term for the UV subdivergence is obtained by applying 
the set of i^-operations identified in the previous step. The individual ii'-operation is carried 
out by the power counting in a similar manner as the ordinal if-operations except that 
the domain of operation is restricted to the corresponding component diagram that the 
subdiagram belongs to. When there are more than one subdiagram, the i^-operations are 
applied successively. The order of applications to different components is indifferent to the 
result. For a proper definition of the integrand Feynman cutoff must be introduced for 
UV-divergent subdiagrams. Of course, the cutoff dependence cancels out in the end. 



3. Building Blocks 

For each component diagram, the building blocks for the subtraction term of the UV 
subdivergence are obtained in the factorized form for the respective UV limit. Recall that 
the building blocks of the IR subtraction term for the annotated forest / is given as those 
obtained for the UV subtraction term associated with the same forest. Thus they reduce to 
the building blocks Bij and U in the UV limit of the forest /'. 



E. Implementation 



We imp 
FORM 



emented the steps described above as separate Perl programs that use internally 
151 ] and Maple. These symbolic manipulation programs take traces, project out 
the magnetic moment, perform analytic integration over momentum variables by means 
of home-made integration tables written in FORM, carry out inversion of matrices which 
creates Bij and U, and execute i^'-operations. The programs for the IR subtraction part are 
integrated with the programs that generate the codes for UV-finite amplitudes previously 
developed [lO|, to form the automated code-generating system for the finite amphtude free 
from both UV and IR divergences. 

The code-generating system takes the name of the diagram and finds the corresponding 
single-line expression of the diagram from the table prepared beforehand. Then it generates 
the numerical integration code in FORTRAN format that consists of the unrenormalized 
amplitude and appropriate subtraction terms. The FORTRAN code is readily integrated 
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FIG. 9: Flow of automated code generation. 



by a numerical integration program, to which purpose we use VEGAS HI], an adaptive 
Monte-Carlo integration routine. 

The whole steps are controlled by make utility and a shell script. It is noted that some 
steps are independent with each other so that they can be performed simultaneously. By 
specifying appropriate options to make, those operations are performed in parallel. This 
feature is beneficial in the present multi-processor/multi-core computer environment. 



VI. CONCLUSION AND DISCUSSION 



In this article we described a systematic procedure for the construction of IR subtraction 
terms in the calculation of lepton anomalous magnetic moment. Our present concern is 
the type of diagrams without internal lepton loops (called q-type diagrams) that account 
for a significant fraction and the most difficult portion of calculation of QED corrections 
to lepton g — 2. We presented a procedure for the automated generation of the numerical 
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integration code, which, together with the previously developed automation scheme for UV- 
finite amplitude, provides prescription for the fully-renormalized and finite amplitude. 

We implemented the code-generating system based on our procedure as a set of Perl 
programs that work with the helps of symbolic manipulation systems, FORM and Maple. 
The entire flow of process is managed by make utility and a shell script. It generates 
numerical integration code in FORTRAN format, which is readily integrated by numerical 
integration system, VEGAS. Our implementation is applicable to any order of diagrams of 
q-type. 

For the purpose of debugging of our automation system we have applied it to the evalu- 
ation of known diagrams of sixth- and eighth-orders. The sixth-order contribution of q-type 
diagrams comes from fifty vertex diagrams, which are reduced to eight self-energy-like di- 
agrams taking account of the Ward-Takahashi identity and the time-reversal symmetry. It 
takes just one minute to create all eight FORTRAN programs for MgQ, (a = a,b, . . . ,h) 
on hp's Alpha machine. Numerical evaluation was carried out on RIKEN's cluster system 
(RSCC). The result obtained is, after incorporating the residual renormalization, 

Af^(q-type) = 0.9052 (11), (49) 

which reproduces the exact result [l6j quite satisfactorily 



4'^(q-type),_, = 0.904979.... (50) 

The value (H9l) was obtained by computation by VEGAS of 2 - 6 hours wall-clock time with 
16 Xeon-CPU's for each diagram, which sampled one hundred million points per iteration 
and 450 iterations. 

The eighth-order q-type diagrams consist of 518 vertex diagrams, which are reduced to 
47 self-energy-like diagrams with the help of the Ward-Takahashi identity and the time- 
reversal symmetry. The entire 47 program sets are generated by our automated code- 
generating system in less than ten minutes on hp's Alpha machine. The numerical evaluation 
is, however, quite non-trivial and requires huge computational resources. For the preliminary 
evaluation we have used 64 to 256 CPU's per diagram for a few months on RSCC to reach 
the relative uncertainty of about 3 %. 

One unexpected byproduct is that it revealed an inconsistency in the treatment of IR 



subtraction terms in the previous calculation of the eighth-order q-type diagrams 
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17|. With 



this inconsistency resolved, the new and old calculations agreed within the numerical preci- 
sion employed. Note that this is the first time that the complete eighth-order contribution 
to g — 2 has been calculated by more than one independent method. The agreement of two 
independent calculations puts the eighth-order contribution to on a firm ground, and 
advances the test of QED and the determination of the fine structure constant derived from 
Oe to a higher level of precision. 

These tests have confirmed the validity of our automation system, and encouraged us to 
tackle the evaluation of the tenth-order contribution with confidence. Of 6354 q-type vertex 
diagrams that contribute to the tenth-order term, 2232 diagrams that have only vertex 
subdiagrams have already been evaluated |10,]. For the remaining 4122 diagrams that have 
self-energy subdiagrams and thus suffer from the infrared singularity, we are evaluating their 
Ward-Takahashi version (254 diagrams) by the automation system described in this article. 
It seems thus far that the numerical results of these diagrams are finite and our subtraction 
scheme works as expected. Some diagrams show somewhat unstable behavior during the 
evaluation which seems to indicate the presence of a severe digit-deficiency problem"^. To 
deal with such cases we have modified the code so that we can carry out the integration 
with extended numerical precision. 

It is important to note that our new subtraction scheme for IR divergences differs from the 
scheme described in Ref. |9| by finite quantity. In the previous approach, the IR subtraction 
term was chosen for a self-energy subdiagram S in the form of Eq. (fT3ll defined by the 
/-operation alone. In the new scheme, we have chosen the form i^^^^) for the subtraction 
term defined as the residual part of L7^(fc) other than the UV divergences. This choice of 
term includes a finite part AL in the previous calculations. Thus the subtraction of the 
nested IR singularity is carried out differently. The latter part of Eq. f|T3|) is related to our 
residual self-mass subtraction, but the form of term Mn* [I] was chosen through the careful 
inspection of the divergent structure of the particular diagram. 

Therefore the previous scheme and the new scheme yield different integrals for diagrams 
containing self-energy subdiagrams. The difference can be traced analytically. In the eighth- 
order case it is confirmed by numerical calculation. The difference is, of course, to be 
compensated by the different set of terms for the residual renormalization step. 



For a detailed discussion of digit-deficiency problem, see Appendix B of Ref. 18|. 
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In this series of articles we have focused on the particular type of diagrams that have 
no lepton loops (q-type). However, our procedure is readily applicable to other types of 
diagrams that are obtained by inserting vacuum polarizations to q-type diagrams. Such sets 



of diagrams of the tenth order are 



io| 



• Set III (a) obtained from the sixth-order q-type diagrams by inserting two second-order 
vacuum polarizations, 

• Set Ill(b) obtained from the sixth-order q-type diagrams by inserting a fourth-order 
vacuum polarization, 

• Set IV obtained from the eighth-order q-type diagrams by inserting a second-order 
vacuum polarization. 

The automated code generation for these sets is accomplished with only small modifications 
to the code. The integration codes have already been obtained, and the numerical evaluation 
has been carried out. The result will be reported elsewhere. 
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APPENDIX A: IDENTIFICATION OF IR SUBTRACTION TERMS: AN EX- 
AMPLE 

As a demonstration of the IR subtraction scheme, we present the identification of IR 
subtraction terms for an eighth-order diagram, Mig, shown in Fig. [TDl 

The diagram Mig has three self-energy subdiagrams denoted by Si, S2, and ^3. The 
number of combinations of those subdiagrams (forests) is seven, and the number of ways to 
assign the types (M, L^, or 6m^) to the nodes is 19. Dropping the cases that contain the 
second-order residual self-mass subtractions 6m^ = 0, it is found that the IR subtraction 
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FIG. 10: An eighth-order diagram Mig. The numerals denote the indices of lepton hnes and 
the roman alphabets denote the indices of photon lines. The diagram Mig has three self-energy 
subdiagrams denoted by Si, S2, and ^3. 

terms are given by the following six cases. On the left-hand side is shown the annotated tree, 
and the figure on the right-hand side expresses the reduced subdiagrams of the component 
terms. 

As a convention we denote the subtraction term by the symbol I and/or R with suffixes 
that refer to the indices of lepton lines. They are the indices to lepton lines contained in 
for /-subtraction, and the range of indices to lepton lines indicated by the left-most and the 
right-most indices contained in 6171^ for /^-subtraction. 
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